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ABSTRACT 

Total and internal energy requirements set constraints on the allowable tem- 
perature on the disk. In a gas pressure dominated disk, they lead to a cooling 
rate that decreases with temperature. These findings are used to show that the 
flow is mildly supersonic, with the azimuthal Mach number moderately exceed- 
ing one in a large extent of the disk. Only in the narrow outermost region, of 
about 12% of the disk extent, the Mach number greatly exceeds one. It is also 
shown that, under most favorable conditions, radiation can , at the most, trans- 
port 20 % of the heat produced, close to the inner radius. Far away, under the 
same conditions, this efficiency can rise to 100 %. This result depends only on 
boundary conditions, being independent of the accretion rate and mass of the 
central object. The geometrical thin disk approximation is violated everywhere in 
the disk, except in the narrow region where, under the unlikely assumption of 
thermodynamical equilibrium or, at least, local equilibrium, the flow is indeed 
supersonic. 

Subject headings: accretion, accretion disks - hydrodynamics, hydrostatic equi- 
librium - Mach Number, supersonic flow - radiation - local thermodynamical 
equilibrium 
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1. Introduction 

In the last 30 years, accretion disks have played a major role in theoretical astrophysics. 
Part of this is due to the growing evidence of their presence in almost every astrophysical 
environment. However, it would be unfair do not recognize that this popularity is also due to 
physical picture of an accretion disk set forward by Shakura & Sunyaev (1973) a-standard 
model. The simplicity of the arguments that have led to the construction of this model 
is quite striking, such that, despite some drawbacks after the long time in the road, some 
authors keep using it as their main tool to explain data. 

However, the model having thinness and hydrostatic equilibrium as one of the main 
assumptions in its foundations, was meant to apply under very specific conditions. These 
certainly would exclude disks with very high accretion rates, or radiative pressure dominated 
disks. As a matter of fact, in that situation the disk would be thermal and secularly unstable 
(Lightman & Eardley 1974). According to the picture given by the standard model, the 
region exposed to swelling or instability would be the innermost one, which would exist 
whenever the radiation pressure exceeds the gas pressure. In other words, for high accretion 
rates. Another difficulty of the model is its impossibility of reproducing the correct transonic 
nature of the flow close to the horizon (Chakrabarti 1990). During quite a long time these 
two issues have been a major concern of the research in the area. Among models aiming 
to develop alternative accretion disk models, it is worth mentioning Paczynski & Wiita 
(1980) model which replaced the usual assumption of local energy balance by a global energy 
requirement, under the assumption of critical radiative flux. Unquestionably, one of the main 
contributions of this work was the pioneering introduction of a pseudo Newtonian potential 
which simulates relativistic effects without any resource to a burdensome mathematics. At 
that time, the development of a theory of geometrically thick accretion disks around black 
holes has been undertaken by Abramowicz et al. (1978); Kozlowzki et al. (1978); Jaroszynski 
et al. (1980); Begelman & Meier (1982). Progress, concerning the transonic nature of the 
flow, came in a much slower way through the contributions from Kozlowzki et al. (1978); 
Liang & Thompson (1980); Paczynski & Bisnovatyi-Kogan (1981); Loska (1982); Muchotrzeb 
(1983); Abramowicz & Kato (1989). Since the flow has to be continuous, the solution that 
incorporates the transonic nature of the flow, close to the horizon, should match the solution 
far away, where, supposedly, the disk is geometrically thin and Keplerian. However, this 
matching introduces some arbitrariness, such as matching radius, and the matching of the 
radial and angular velocities. The problem here is the angular momentum distribution, 
which, close to the horizon, is highly dependent upon the matching radius (Chakrabarti 
1996). Once this radius is determined, the angular momentum follows practically constant 
inside the horizon, a fact which suggests null boundary condition assumption for the torque 
(Muchotrzeb 1983). Besides, if one adopts a more realistic model of disk, allowing for 
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dissipation, the situation may become much more complex. This is due to the topology of the 
solutions near the critical and subcritical points, both of them occurring in dissipative flows. 
Critical points satisfy regularity conditions in the form of algebraic constraints, while for 
the subcritical points they are given by integral requirements. In determining the character 
of the solution near these points, the role of the viscosity is quite decisive, such that there 
is a critical value of a above which no stationary accretion is possible (Muchotrzeb 1983), 
a question contested by Matsumoto et al. (1984) and Abramowicz & Kato (1989). Clearly, 
the matching of an "inner disk" with an "outer disk" is highly dependent on what kind of 
disk one assumes as an outer disk, that is to say, geometrically thin or not, a fact that may 
be seen in the works of Liang & Thompson (1980) and Loska (1982). 

Recently, these issues have been revived, connected with the question of "no torque" 
inner boundary condition (Krolik 1999; Gammie 1999; Agol & Krolik 2000; Afshordi & 
Paczynski 2003). As argued by Afshordi & Paczynski (2003), the no torque inner boundary 
condition applies whenever the disk thickness is small everywhere. However, results based on 
semianalytical and numerical work employing advective dominated accretion flows (ADAF) 
and convective dominated accretion flows (CDAF) show that accretion disks are thick for 
very high and very low accretion rates (Ball et al. 2001). Thinness would occur only for 
moderate accretion rates (0.01 < L/L Edd < 0.1). It seems that, in order to bypass the issue 
of thermal balance, one adopts the simplifying assumption of disk thinness. This is not 
justified a priori, because even numerical simulations yield relatively thick outer boundary, 
not far away from the inner one (Afshordi & Paczynski 2003). However, these simulations 
suffer from limitations, one of them being not including cooling processes. In the literature, 
alledgedly, the assumption of highly supersonic character of the Keplerian flow is used as the 
argument for the disk thinness. This is, for sure, intuitive, but, by no means, unquestionable. 
Somehow, it seems that thinness is favored by the usual assumption of local equality between 
heat production and cooling by radiation, which overlooks the problem of energy conservation 
in the disk. One of the advantages of doing so is to render the problem algebraic, rather 
than differential , non local. In the solutions to the disk equations, the a parameter doesn't 
enter with a high power (Frank et al. 1992), which due to its inherent uncertainty, is very 
interesting. At this point, we realize how decisive the use of the geometrical thinness, as a 
boundary condition, is in determining the structure of the whole disk. If cooling of the disk 
is not efficient, most of the viscously produced heat is stored as proton internal energy and 
advected with the flow. The temperature tends to the virial value and the height scale of 
the disk will be comparable to the radial distance. More important, the internal energy of 
the plasma approaches the rest energy and, at the inner radius, the sound velocity is much 
greater than the radial velocity, being a substantial fraction of the light velocity. Since the 
gas only cools after passing through the sonic point (Muchotrzeb 1986), the sonic velocity 
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keeps growing till the sonic point which, in that situation, will be closer to the hole. Though 
the picture we have given may be incomplete, it seems that, somehow, the importance of the 
disk being thin or not has been overlooked. To justify this very pertinent point, we would 
like to stick on the conclusions by Afshordi & Paczynski (2003) and on the possibility, cast 
by them, of non existence of thin geometrical disks in nature. We, therefore, propose to 
study, in a more detailed way, the geometrical thin disk assumption. Our procedure will 
be to use, from the very beginning, energy requirements constraining both the temperature, 
the height scale of the disk, viscous heating, and cooling of the disk. This, certainly, is a 
global requirement that avoids the usual assumption of local equality between viscous heat 
production and cooling. Since the picture of geometrical thin disk is so profoundly rooted 
in our minds, in this work we shall focus our attention only on that. 

The plane of the paper is as follows. Section 2 presents the system we are considering, 
the assumptions and the disk equations, including an outline of the ingredients to relate 
angular momentum transport and dissipation to the density and radial velocity. In section 
3, the total and internal energy equations are handled in such a way as to lead to an equation 
to the temperature of the disk. In section 4, some approximations are made in order to solve 
the equation for the temperature in conditions similar to those where the standard model 
is supposed to apply; the main differences between our solution and the standard one are 
emphasized. In section 5, general considerations about the radiative transfer problem link the 
results to specific radiative mechanisms. In section 6, the flow regimes that may emerge are 
discussed in the context of different boundary conditions. In section 7, the main conclusions 
are drawn. 



2. The System And The Disk Equations 

The system we consider is a disk around a primary compact object (a neutron star, a 
black hole), of mass M. Self gravitation in the disk is negligible. Working in a reference 
frame tied to the center of mass of the compact object, a point on the disk is at a distance 
R from the center and has coordinates (x, y, z) or (r, 0, z). Of course, R = \J (r 2 + z 2 ). 

For the moment, we shall use Newtonian gravitation, and the potential at the point (r, 

0, z) is 



_ GM 

G being the gravitational constant. 
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For the <p component of the velocity we will assume a Keplerian profile. Therefore, we 
shall write 



V <j) = rQ k , 

where fl k is the Keplerian angular velocity, 



Since we are assuming the disk in hydrostatic equilibrium, 

dP 

dz 

where p is the density and P is the pressure. Now, writing 



P = P g + P r , 

where 

P 9 = P — T 

m H 

is the gas pressure, T is the temperature, and 

F r p a T £ 



P 



c ran 

is the radiative pressure. F r is the radiative flux, £ is the semi-height scale of the disk, k B 
is the Boltzmann constant, m H is the Hydrogen mass, c is the velocity of light, and a T is 
the Thompson cross section for electron scattering. Assuming constant density along z, and 
taking the average over z in equation (1), yields 

n k 2 f _ k B T | F r a T e 

3 mi c mi 
The continuity equation is written as 



M = -4 7T r £ V r p , 



(3) 
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where M is the accretion rate, and V r is the radial velocity. 

Now, we need to work out some relations between angular momentum transport and 
dissipation to obtain the usual expressions for the density and radial velocity. If J in , J and 
J ou t are, respectively, the inwards angular momentum transport (per unit time), the rate at 
which the angular momentum flows into the central compact object and the outwards flux 
of angular momentum, we have 



and we may write 



Jin — Jo + Jouti (4) 

j in = M y/GMr, (5) 
Jo = M y/GMn (6) 



and 



j out = 2irr £r r(f) , (7) 
where r r ^ is the stress tensor given by 

T r<t> = -2rja r}(j) , (8) 
a r ^ being the rate of strain tensor, given by 

r dVL k 

°r>+ = -2-fr- (9) 

Then, the angular momentum conservation may be written as 

M VG Mr = M y/GMn + 2Txr ir r<S) . (10) 

In the above equation, we have assumed null boundary condition for the torque at the 
inner radius r = r\. This yields for the energy dissipation, 



-7- 



and for the density 



» U-r^T). (12, 



8tt n k P V V r 



3. Energy Conservation In The Disk 

The equation for the total energy conservation in the disk, as the matter flows from r 2 
to r, may be written as 

P \ \ v i + i] v k- PSikVt] = -pv k - A, (13) 



dxk V \2 / / dxk 

where £ is the total internal energy, A is the amount of energy per unit volume that leaves 
the system. Assuming the system is cooled by radiation, we may write 

A = £■ < 14 > 

Inserting this expression into equation (13) and using the equation of continuity, we 
obtain 

M Q^ + e + ^T) =M^L-Atc J\F r dr. (15) 

Finally, assuming uniform viscous heat production along z, we write the internal energy 
equation as 

d ( C\ D dVk dFr , 1 + ftR\ 

dt ( ^ ) = " p ^- dT + T- (16) 

Now, making the approximation 

~d~z~ ~ T' 

we solve for F r in equations (15) and (16), and equate the results to obtain 

d , Pi d , ^ M d /1_ 2 tt k 2 l 2 \ M GM 



-eV r -(pO (rV r )+q + = -— - -ltf + £ + ^L_ _ (17) 

ar r ar 4 7r r Or V 2 6 / 47r r 
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This is the final energy equation for the disk, which has to be solved, together with the 
hydrostatic equilibrium equation, to obtain the temperature profile along the disk. It should 
be remarked that its deduction has been achieved with no further assumption concerning 
the flow regime, i.e., supersonic or subsonic. We have, however, neglected terms of order of 
V 2 , and haven't considered transport of radiation along the radial direction. 



4. The Solution To The Energy Equation 

In the following we shall be solving equation (16) with the main concern to find out 
if there are geometrical thin solutions for the disk. It is believed, under the a-standard 
prescription, that far away from the inner edge, in region dominated by gas pressure, the 
disk is geometrically thin. Therefore, the following analysis is meant to this region. In other 
words, we shall assume r >> r±, and P g » P r , which implies 

2/>2 



and 



k B 

~ 1 , 

3 m H 



n 2 e 2 



Therefore, we may write 



Pi 



3m H \ 1/2 MQ k 2 



2 6 tt 2 AW aTV2 ' 



V r ~ a 



and 



m H r Q k ' 

Pi w -^—MQ k , 
An a 

( h \ 3 ^ 2 T 3 / 2 
\m H J rVL k 



, k B \ T 
r V r ~ — a 



m H J Q k ' 



which, inserted into equation (17), yields 

dT _UT _8m H GM 
dr 11 r llks r 2 ' 



(18) 
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the solution being 

_ 14 m h G M . . 

T = Arn +0.32— . 19 

k B r 

Using equations (15) and (19), we may write for the radiative flux, 

_ 3 GMM _ 70 k B M 

r 40 vr r 3 d>d>itm H r*/ u 1 ' 

Equation (20) can be cast in a much more interesting form if we use equation (19) to 
express the constant A. We, then, obtain 

_M_ ( 29 GM 70 k B \ 
r 8nr 2 \11 r llm H J ' 

or 



M /29GM_70 y 
8nr 2 \11 r 33 k ) K J 

The first term in this equation is, practically, the same one that comes from usual 
accretion disk theory, i.e., 

3 GMM 



8 7r r 3 

The second, that doesn't appear in the usual theory, displays the effect of the temperature, 
or the height scale: the hotter, or the thicker, the less it emits. It should be remarked that 
the largest temperature is 

29 m H GM 



Tm — 



70 k f 



In the standard formulation, one assumes equality between heat production and cooling 
by radiation. Since the energy produced due to viscous dissipation, q + , is independent of the 
height scale, no matter how thin or thick the disk may be, it always adjust itself to get rid 
of this energy. As compared to this, the present formulation constitutes an improvement: it 
is clearly seen that, if the disk thickens, it starts dimming. 
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5. Radiation, Electronic Temperature and Height Scale 

If we write the energy exchange between protons and electrons as 

3 pk (T\ - T e ) 
ep ~2m H t ep ' 

and equate it to the flux given by equation (21), we obtain the time protons and electrons 
need to reach thermal equilibrium, 



tep ~ ^ ( 1 rp ) tdyn i 



tdyn being the dynamical time. We shall assume a « 1, and t ep » td yn - 

The temperature we have obtained, due to our starting assumption, has to be identified 
with the ion temperature. Clearly, this conclusion breaks down when there is an abundant 
pair production and the pairs dominate the pressure. In such a situation, it has to be 
identified with, what we improperly call, the equivalent temperature, which we shall define 
at the right time. To know the electronic temperature, we should assume a specific radiation 
mechanism. Therefore, in the following, we shall make general considerations to link our 
previous results to the radiative transfer problem. In our approach, we shall try, to make a 
formulation not depending, as much as possible, on the values of the accretion parameters. 

Our starting point is to know the optical depth for thermal bremsstrahlung absorption, 
r ab . By definition (Frank et al. 1992), 

r ab = 6.22 x 10 22 p 2 T- 7/2 £. 

We should expect r a b decreasing as we come close to the hole. Then, in the inner region, the 
disk emits thin bremsstrahlung, and the flux is (Rybicki & Lightman 1979), 

F = 5 x lOVT' 1 /^. (23) 

This yields 

M 9 

r- « 1.38 X lO-^r 5 

(r in units of ^-). 

Finally, we obtain the electron scattering optical depth, 

0.225 Mi 7 _ 1/2 
a Mio 
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We conclude that if 

Mn < 3M 10 , 

in such a way as to have the gas dominating the pressure, there may be three regions in the 
disk, depending on the accretion parameters, a, M, M: 

a) The innermost region: r sc > r a b, and r a b < 1. In that region, although T e < Tj, it 
may be large enough to produce pairs. If the pairs are very numerous, they may dominate 
the pressure. Then, we should replace Tj for the equivalent temperature, i.e., 



T eq = T i (l + {l + 2y) 



where y is the ratio pair density over ion density, that implies, 

^fc^ 2 = T eg . 
m h 

To obtain the electronic temperature, we have to correct equation (23) for the inclusion 
of pairs, i.e., 

F = 5 x 10 20 p 2 T- 1/2 £ (1 + 2y) . 
This will lead to an equation relating T e and i, i.e., 

(IT "33 \R' ' ' (2J) 



where 

M 2 

A = 5.13xl0 41 - i?F (l + 2 2/ ) 2 T e . (25) 



Clearly, for a given T e , we may interpret it as an equation for £, which will have solution 
as long as A is less than a critical value A c . Solving under critical conditions, we obtain, 

6.95 - 11.185a; 2 + 3.78a; 4 = 0, (26) 

with x — -|, the solutions being x = 0.942 and x = 1.439. The solution x = 1.439 leads to 
negative flux. Therefore, only x = 0.942 is acceptable. 

Equating A, given by equation (25), to A c , given below, 

^c=(|) 2 ^ 1 V 4 , 
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we obtain T e . Therefore, in that region 

^ 0.942i?. 

b) The intermediate region: r sc > r a b > 1. In that case, the disk will emit modified 
black-body-like, and the flux is (Rybicki & Lightman 1979) 

F MB = 2.3 x 10y/ 2 T e 9 / 4 . (27) 
Now the equation relating T e and £ will be 



A _ / 29 70 
£3 ~ I IT ~~ 33 Vi? 



2 x 2 



with 



A = 7.98 x 10 : 



9/2 



2 



Proceeding in the same way, as we did before, we obtain 

6.95- 11.185a; 2 - 3.496a; 4 = 0. (28) 
The solution now is x = 0.730, the critical A is 



/280\ 2 . 

A c = R 3 x 7 , 

V 99 J 

and 

0.730i?. 

c) The outer region: r a 5 > 1 and r a b > r sc . We shall assume that Tj = T e . If we assume 
the disk emits black-body-like, then 

-- + -x 2 + Ax s = 0, (29) 

where 

A » 6.1 x 10 15 Ml2 ° 



M i7 r 

For r = 1000, M 10 = M i7 = 1, 

x w 0.03. 



It should be remarked that in this region, there is no critical A, and 

0.03.R. 
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6. Boundary Conditions and Flow Regime 

The expression for the radiative flux may be rewritten, 



MK 2 /29 , 70' 

m 4>-^) » ( 30 ) 



8vr V 11 33/ ' 

where V s and are, respectively, sound velocity and azimuthal Mach number. Of course, 

> 0.9. 

Our problem, now, is to know how much the azimuthal Mach number can exceed the 
above value. For that, let us assume that at the outer radius, R d , the temperature is T d . 
Then, 

T,- 0.32^^1) (31) 
If ( is the efficiency for converting gravitational energy into radiation, 

(Mc 2 = 0.3^-0.4^ , (32) 

n m H \ k B nr d J 

or 

£("-")-<)■ 

ri, and are, respectively, the inner and the outer radii in units of half gravitational radius, 

Rg_ 

2 • 

Using now the definition of the Mach number and the expression for the temperature, 
we obtain 

o 3 kn R^ 

m J = l-_«_ —A, 34 

which allows us to conclude that if > 1, A < 0. Since, T d > 0, we must have ( < If 
( is maximum, 



(35) 
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which, inserted into equation (34), leads to 

Therefore, for maximum radiation efficiency, there will be a region in which the flow 
regime is supersonic. However, even in that case, the disk will be moderately supersonic in 
most of its extent. Just in a small outer fraction, which depends on the size of the disk, will 
be really supersonic. It should be remarked that maximum efficiency means null temperature 
(or almost) in the outer border of the disk. To make things clear, let us calculate for a disk 
with r d = 200, for: (a) r = 195 and (b) r = 199. We obtain: 



r = 195 -> = 4.23 
r = 199 -> = 9.40 

For different boundary conditions, such that A > 0, the flow regime will be, practically, 
sonic in the whole extent of the disk. 

We have, here, a situation completely different from that we encounter in the standard 
model. There, due to the equality between dissipation and cooling, radiative flux obeys the 
same boundary condition imposed on the torque, and has a maximum at R — 1.36 R\. 

Inserting the value for A, obtained from the luminosity constraint, in the expression for 
T, given by equation (19), we obtain, at r = r±, 



R 



0.236 



- = 0.106. 
R 



T ~ 5.84 x 10 11 , 

independent of the mass and accretion rate. For /3 ~ 1, this implies, for the azimuthal 
velocity, a Mach number of about 1. Therefore, close to the inner radius, 




With these values, 

p 

-k5x 10~ 18 M, 
p 

and our estimates hold as long as M < 2 x 10 17 gs~ 1 , for M 10 = 1. It is interesting to remark 
that the ratio of the radiative flux to the heat flux due to viscous dissipation, close to the 
inner radius, is 
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F r 



1 



5 



It should be remarked that, if some other cooling mechanism is considered (Meirelles Filho 
et al. 1984), this could be even less. 

Analyzing the r dependence of the radiative flux, we see that, asymptotically, under the 
same boundary conditions, 




the importance of radiative cooling grows as we go towards the outer radius. Therefore, 
the geometrically thin disk assumption only holds in a very narrow region where the flow is 
supersonic. As a matter of fact, the results we have just obtained invalidates the geometrical 
thin disk assumption. 

We should insist in a point: this region, where the flow is indeed supersonic, only exists 
if radiation efficiency is maximum. In our formulation, ( = 0.05 and ( = 0.3, respectively, 
for a Schwarzschild and a Kerr hole. In that case, 

r < 0.88r d -> m < 2.0, 
r > 0.88r d -> m > 2.0. 

If it is a little less, A > and the flow will be sonic, practically, all over the disk. 

Finally, it should be said that, since the radiative flux continuously grows as we approach 
the hole, it can not satisfy the same boundary condition the torque does. Therefore, no 
matter the correct boundary condition, the hole does swallow, always, an amount of energy 
per unit time, which is, at least, of the order of the disk luminosity. 



We have used total and internal energy equations to obtain an equation for the tem- 
perature along the disk. This equation solved under gas pressure dominance shows a disk 
hotter than the standard one. Besides, it implies a radiative flux decreasing with tempera- 
ture or height scale. Though not considered here, these findings impose severe constraints 
on the development of instabilities in the disk. In that formulation, the role of advection is 
more important than usually admitted. Independently of mass, accretion rate, and viscosity, 
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Conclusions 
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radiation transports, at the most, 20 % of the produced heat, close to the inner radius. Far 
away, its efficiency may rise to 100 % . Under the assumption of large radiation efficiency, 
as a boundary condition on the solution for the temperature, it is shown that the flow is 
supersonic, the azimuthal Mach number moderately exceeding lin most of the extent of the 
disk. Only, on a narrow outer region, of about 12 % of the disk extent, the Mach number 
greatly exceeds 1. This highly supersonic regime is only obtained under the assumption of 
TE or, at least, LTE, which is very unlikely to occur (Hubeny et al. 2000, and references 
therein). These results show that the geometrical thin disk assumption fails everywhere in 
the disk, but in the outer narrow region. However, even if the assumption of LTE holds in 
that region, it represents a very small part of the disk. 

REFERENCES 

Abramowicz, M., Jaroszynski, M., & Sikora, M. 1978, A&A, 63, 221 
Abramowicz, M., & Zurek, W.H. 1981, ApJ, 246, 314 
Abramowicz, M. 1985, PASJ, 37, 727 

Abramowicz, M., Lanza, A., & Percival, M.J. 1987, ApJ, 479, 179 

Abramowicz, M., Czerny, B., Lasota, J.P, & Szuszkiewicz, E. 1988, ApJ, 332, 646 

Abramowicz, M., & Kato, S. 1989, ApJ, 336, 304 

Afshordi, N., & Paczynski, B. 2003, ApJ, 592, 354 

Agol, E., & Krolik, J.H. 2000, ApJ, 528, 161 

Ball, G.H., Narayan, R., & Quataert, E. 2001, ApJ, 552, 221 

Begelman, M.C., & Meier, D.L. 1982, ApJ, 253, 873B 

Chakrabarti, S.K. 1990, Theory of Transonic Astrophysical Flows, World Scientific 
Chakrabarti, S.K. 1996, Phys. Rep., 266, 229 
Czerny, B. 1986, Acta Astronomica, 36, 1 

Franck, J., King, A., & Raine, D. 1992, Accretion Power in Astrophysics, Cambridge Uni- 
versity Press 

Gammie, C. 1999, ApJ, 522, L57 



-17- 

Hubeny, L, Agol, E., Blaes, O., & Krolik, J.H. 2000, ApJ, 533, 710 

Jaroszynski, M., Abramowicz, M., & Paczynski, B. 1980, Acta Astronomica, 30, 1 

Kozlowzki, M., Jaroszynski, M., & Abramowicz, M. 1978, A&A, 63, 29 

Krolik, J.H. 1999, ApJ, 515, L73 

Liang, E.P., & Thompson, K.A. 1980, ApJ, 240, 271 

Lightman, A.P., & Eardley, D.M. 1974, ApJ, 187, LI 

Loska, Z., B. 1982, Acta Astronomica, 32.1.2, 14 

Matsumoto, R., Kato, S., Fukue, J., & Okazaki, A.T. 1984, PASJ, 36, 71 
Meirelles Filho, C, Lima, C.L., Miyake, FL, & Timoteo, V. 2003, ApJ, 595, 322 
Muchotrzeb, B., & Paczynski, B. 1982, Acta Astronomica, 32.1, 1 
Muchotrzeb, B., 1983, Acta Astronomica, 33, 79 
Muchotrzeb, B., 1986, Acta Astronomica, 36,1 
Paczynski, B., & Wiita, P.J. 1980, ApJ, 88, 23 

Paczynski, B., & Bisnovatyi-Kogan, G. 1981, Acta Astronomica, 31.3, 283 

Rybicki, G.B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics, John Wiley and 
Sons, New York 

Shakura, N.I., & Sunyaev, R.A. 1973, A&A, 24, 337 



This preprint was prepared with the AAS IATgX macros v5.2. 



